Coarse grained models of stripe forming systems: phase diagrams, anomalies and scaling hypothesis 
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Two coarse-grained models which capture some universal characteristics of stripe forming systems are stud- 
ied. At high temperatures, the structure factors of both models attain their maxima on a circle in reciprocal 
space, as a consequence of generic isotropic competing interactions. Although this is known to lead to some 
universal properties, we show that the phase diagrams have important differences, which are a consequence of 
the particular k dependence of the fluctuation spectrum in each model. The phase diagrams are computed in a 
mean field approximation and also after inclusion of small fluctuations, which are shown to modify drastically 
the mean field behavior. Observables like the modulation length and magnetization profiles are computed for the 
whole temperature range accessible to both models and some important differences in behavior are observed. A 
stripe compression modulus is computed, showing an anomalous behavior with temperature as recently reported 
in related models. Also, a recently proposed scaling hypothesis for modulated systems is tested and found to be 
valid for both models studied. 
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I. INTRODUCTION 

Stripe forming systems are common in nature. Examples 
range from solid state systems, like ultrathin ferromagnetic 
films^, stress-induced submonolayer islands on several sub- 
strates^ and strongly correlated electron systems^, to soft 
matter systems like Langmuir monolayers^, block copoly- 
mers^, colloids and soft core systems^ili. The origin of 
stripe structures can be traced back to some strong univer- 
sal mechanisms^. At the heart of stripe formation there is 
usually a competition between different kinds of interactions 
which lead to frustration and the formation of patterns or mi- 
crophases. Common examples involve competition between a 
short range attractive interaction and a long range repulsion-^ 3 .. 

Although several general properties of stripe forming sys- 
tems are well known, as mean field phase diagrams and the 
behavior of the order parameter near the transition tempera- 
ture^—, the behavior at lower temperatures is less explored 
and new phenomena continue to appear. An example of this is 
the recently reported reentrance in the phase diagram of frus- 
trated dipolar ferromagnets in an external magnetic field 17 ' 18 . 

These facts motivate us to look for a deeper understanding 
of the phase behavior of these kinds of systems. In the present 
work we address several questions which to a large extent re- 
main open. We present a complete study of two very well 
known models of stripe forming systems, defined at a coarse- 
grain level, in which the main characteristic defining the phase 
behavior of the systems is the form of the effective spectrum 
of fluctuations, i.e. the wave vector dependence of the bare in- 
verse correlation at low energies. This is known to determine 
the existence and type of phase transition, at least at mean 
field level. We define a coarse-grained free energy which is 
suitable for analyzing the low temperature behavior of rele- 
vant quantities, like magnetizations and modulation lengths, 
besides allowing the determination of the complete phase dia- 
grams, in the whole temperature range. We show results from 
mean field approximation and also when fluctuations are taken 
into account. An efficient numerical implementation of the 



solutions allows us to compute several quantities with high 
accuracy. 

The two models studied here are usually considered in the 
literature as yielding essentially the same physics. Instead, 
we show that besides obvious common characteristics, both 
models have very different phase diagrams and very different 
behaviors of important quantities, like the temperature depen- 
dence of modulation length. By introducing fluctuations to 
the mean field results we show that the nature of the transitions 
changes and the behavior of key quantities like magnetizations 
are drastically affected by fluctuations. An important charac- 
teristic of stripe forming systems with long range interactions 
is the temperature dependence of the modulation length and of 
the domain wall widths. We compute the T dependence of the 
modulation length in the whole T range, verify the well known 
square dependence near the transition and show that our ana- 
lytic results compare very well with recent experimental data 3 
down to low temperatures. We also introduce a measure of 
the width of domain walls and compare the T dependence of 
the width relative the behavior near the transitions, where a 
single mode approximation usually works fine. Finally, mo- 
tivated by recent work on a scaling hypotesis for modulated 
systems^ we compute a response function, the stripe com- 
pressibility, and show that it presents an anomalous behav- 
ior at low temperatures, in agreement with scaling predictions 
of systems with power law competing interactions. Interest- 
ingly, the general scaling behavior seems to work well also 
for a model which does not correspond to a system with long 
range interactions, our model 2 defined below. This points 
to a wider applicability of the proposed scaling hypotesis for 
stripe forming systems. 

The paper is organized as follows: in Sec. II we introduce 
the coarse-grain free energy in terms of a generic spectrum 
of fluctuations. We define the equations of state at the mean 
field level and when fluctuations are included and show how 
to solve them in an efficient way. In Sec. Ill we introduce the 
two particular models studied, and present the main results 
on the phase diagrams and temperature dependence of several 
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quantities. In Sec. IV we present the results of the stripe 
response function for both models. Finally, the conclusions of 
our work are discussed in Sec. V . 



II. MEAN FIELD APPROACH INCLUDING 
FLUCTUATIONS 



We consider a generic model in two spatial dimensions de- 
fined by the coarse-grained free energy: 
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Here, a represents the lattice spacing, f3 = l/fc^T is the in- 
verse temperature, 7 and 6 are phenomenological constants 
and 0o is the saturation value of the order parameter. The first 
term favors homogeneous configurations of the order parame- 
ter and the second term is intended to represent a competing 
force, represented by the kernel J(f). These two terms can 
be considered the continuum limit of more microscopic com- 
peting interactions, as e.g. exchange and dipolar interactions 
in thin magnetic films^., hydrophobic and hydrophilic com- 
ponents in microemulsions or block copolymers^, or charge 
separation due to long range Coulomb forces in low dimen- 
sional electron systems 5 . The third term represents the mean 
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field entropy. This effective Hamiltonian can be formally ob- 
tained, e.g. by a Hubbard- Stratonovich transformation^ plus 
a coarse-grain process from an Ising model with competing 
short and long range interactions (represented by the non lo- 
cal term). A more common approach assumes an expansion of 
the entropy term for small 0, which leads to the well known 
double well potentia l 15 ' 20 . We will consider instead the full 
form given by ((T), which allows to compute the phase dia- 
grams down to low temperatures. 

Writing the quadratic part of ([!} in reciprocal space and 
normalizing the order parameter field to its saturation value 
yields: 
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where the spectrum of fluctuations is given by 
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A{k)= 1 k 2 + -J{k). 



(3) 



We consider systems in which the fluctuation spectrum has 
a single isotropic minimum at a non zero wave vector fco. 
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If .A (fco) < the system will develop modulated structures. 
With some additional transformations it is possible to express 
the effective Hamiltonian in dimensionless form: 



= / d 2 fci(fc)0£0r + 4 / d 2 x{{l + 0(f)) ln(l + 0(f)) + (1- 0(f)) ln(l- 0(f)) -21n2} 
E Q J fj , 
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In this way all magnitudes but the characteristic energy Eq in 
expression (0]i are dimensionless . This expression was used 
for all the numerical computations in this work. 
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where the average is relative to the mean field measure. 



A. Stripes solutions 

It is well known that, in the absence of external fields, the 
solutions which minimize the free energy of model are one 
dimensional modulations of the order parameter, i.e. stripe 
solutions in d=2, which can be written as: 



4>{x) = TOj sin ^27r(2i 
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(8) 



where A is the "modulation length". In all calculations, the 
number of modes considered was i max = 60. This allow us 
to make a full characterization of the modulation profiles, as 
shown in the following sections. 
The mean field equation (O reads: 



(x) = tanh 



d 2 k 



(9) 



To obtain the solution at a given temperature, equations (|8) 
and (O are solved for fixed A. Then the free energy is mini- 
mized with respect to A to get the final solution. 



B. Stripes solutions including fluctuations 

Up to second order in the fluctuation fields, equation (|7) is 
given by: 
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In order to solve this equation the mean squared local fluctu- 
ations (ip 2 )(x) have to be computed. To do this we used the 
following method: let's consider the field <fi(x) in the pres- 
ence of a local unknown molecular field due to the rest of the 
system. 

The partition function for a single site in an external field h 

is: 
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FIG. 1 . (Color on-line) Local relation between the fluctuations am- 
plitude and the mean order parameter. 
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Of course, the value of the local molecular field h remains un- 
known. Numerically, we can consider it as a free parameter 
to determine parametrically the non trivial relation between 
the mean square fluctuations (i(j 2 }((j)) and the mean local or- 
der parameter. The parametric solution of the system (Q~2} is 
shown in FigQ] It is important to note that, as a result of the 
temperature dependence of the potential used, the relation be- 
tween (ip 2 ) and (4>) is temperature independent. 

The equation of state including the fluctuation fields now 
reads: 
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The effect of fluctuations is to add a new term in the effec- 
tive field, which tends to rise up the order parameter profile 
in the region corresponding to domain walls. This means that 
when fluctuation effects are present profiles should be more 
sharp and correspondingly domain walls will be thinner. 

In the following we apply the method described above to 
solve equations (O and ( fT3l for two well known models with 
stripe phases. 
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III. PHASE DIAGRAMS. 

We focus our study in systems which present isotropic in- 
teractions, meaning that A(k) = A(k). Due to our selection 
of dimensionless variables, the minimum of the fluctuation 
spectrum is reached at k = 1 and then A(l) = —1 . Be- 
sides these universal characteristics, different models are de- 
fined by the particular form of the fluctuation spectrum. The 
first model considered (model 1) is defined b y 14 ' 18,21 : 



A(k) = -l + a(k - l) 2 
and the second (model 2) by 2 ^: 

A(k) = -l + a(fc 2 - l) 2 . 



(14) 



(15) 



Model 1 has a leading k dependence, for small k, which is 
linear in k = \k\, i.e. it is non-analytic. It is a good model of, 
e.g. an ultrathin ferromagnetic film with strong perpendicular 
anisotropy, in which short range ferromagnetic exchange in- 
teraction competes with the long range dipolar interaction — . 
Model 2 is known as the Swift-Hohemberg model, and was 
introduced to describe the physics at a convective instabil- 
ity ^^ _ Alternatively, it can represent the continuum limit of a 
system with isotropic attractive nearest neighbors interactions 
plus repulsive next-nearest neighbor interactions 

Although both models are similar in that they have a single 
minimum at a non zero wave vector modulus, they neverthe- 
less show important differences which are shown in what fol- 
lows 2 ^. To our knowledge a detailed comparison of models of 
this kind has not been done up to now, and this frequently led 
to erroneous conclusions about their behavior. 



A. Model 1 

Fig. |2]shows the mean field phase diagram for model 1 . The 
parameter a is proportional to the curvature of the fluctuation 
spectrum at the minimum wave vector: 



1 
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Note that, due to our choice of dimensionless variables, the 
temperature dependence is included in (3, defined as the in- 
verse of T/T c . This is the reason why the critical line is hor- 
izontal. This line defines a second order phase transition be- 
tween a low temperature modulated phase and a high temper- 
ature disordered one. 

The dashed-dotted curves represent lines of constant mod- 
ulation length A. In fact, for fixed a the modulation length de- 
creases steadily as the temperature increases towards T c and 
then curves nearer the critical line correspond to smaller val- 
ues of A. The full temperature dependence of the modulation 
length is shown in Fig. [4] (dots), which shows that it varies 
continuously with temperature with an asymptotic value cor- 
responding to the minimum of the spectrum. 

The phase diagram of model 1 including fluctuations is 
shown in Fig. [3] At first sight it seems similar to the mean 




FIG. 2. (Color on-line) Mean field phase diagram for model 1 . The 
dashed-dotted lines represent equal modulation length (A) curves . 
The upper full line (blue on-line) corresponds to a second order tran- 
sition between a modulated and a disordered phase. Dots represent 
the results of the numerical solution. 
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FIG. 3. (Color on-line) Mean field phase diagram for model 1 includ- 
ing fluctuations. The dashed- dotted lines represent equal modulation 
length (A) curves . The upper full line (blue on-line) corresponds 
to a discontinuous transition between a modulated and a disordered 
phase. Dots represent the results of the numerical solution. 



field diagram. The critical temperature is depressed relative 
to the mean field one, as expected when fluctuations are in- 
cluded. The main difference with the mean field results is that 
the transition is discontinuous, i.e. fluctuations change the na- 
ture of the phase transition. In this case, the critical line was 
determined at the temperatures where the difference between 
the free energies of the disordered and modulated solutions 
changes sign. This characterizes a first order transition in- 
duced by fluctuations. 

The behavior of the lines of constant modulation length is 
different from that of the mean field case. In Fig. |2]all the lines 
of constant A approach a limiting value of A = 1 at T c , the 
value corresponding to the minimum of the spectrum A(k), 
as can be seen in Fig. [4] (upper dots). When fluctuations are 
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FIG. 4. (Color on-line) Temperature dependence of modulation 
length (A) of model 1 for a = 1, normalized to the value correspond- 
ing to the minimum of the spectrum, Ao. Superior and inferior dots 
correspond to the results of the case without fluctuations and when 
fluctuations are included, respectively. The full lines represent the 
best fits of numerical results by the proposed functionality d!7l l. The 
lower full line (red on-line) is interrupted because of the transition to 
the disordered phase. 



included the transition line is not an equal modulation length 
curve, in such a way that the modulation length at the transi- 
tion point is a decreasing function of a. 

While it is well known that the modulation length grows 
quadratically with temperature near the transitio n 18 ' 27 , the low 
temperature behavior of A is less known. Based on scaling 
approaches from other authors, Mentes et al.- proposed a nu- 
merical fit for their data on stripe forming ultrathin Pd/W(l 10) 
films. The proposed functionality, which takes into account 
fluctuation effects on the interface free energy and domain 
wall width, is the following: 



A(T) = Ao + a(l - bT) exp(-cT). 



(17) 



In Fig. [4] we show the results obtained directly from our nu- 
merical solutions (dots) together with the best fits from equa- 
tion < fT7T > (solid lines) for both approximations, mean field and 
mean field plus fluctuations, showing that the physically mo- 
tivated form ( fTTI i is compatible with our numerical data. 

In figure the behavior of the amplitude of the modulated 
solution is shown. The temperature in the horizontal axis is 
scaled relative to the critical temperature of the mean field so- 
lution. A notable fact is that the first order transition induced 
by fluctuations is very strong, the amplitude stays at very high 
values in the whole modulated phase, until the transition takes 
place. This is also reflected in the behavior of the domain 
wall width relative to the modulation length, which is shown 
in figure [6] The definition used for the domain wall width was 



Jo 4 f( x ) x dx 

So /(*) dx 



(18) 



where the weight function f(x) is f(x) = M — <f>(x), and 
M corresponds to the maximum value of the <f>{x) profile. In 



FIG. 5. (Color on-line) Temperature dependence of the modulation 
amplitude (AT) for model 1 with a = 1. The lower full line (blue on- 
line) corresponds to the mean field solution and the upper full one 
(red on-line) is the solution when fluctuations are included . Dots 
represent the results of the numerical solution. 



our definition of 4>(x) we fix the phase of the modulation, tak- 
ing 0(0) = and <f>(x) positive in the interval (0, A). For 
example, for a sine like profile, the domain wall width is: 



£ 

A 



/ o *(l-sin(0))fldfl 
7r/*(l-sin(0))d0 



0.1303 



(19) 



This value would be a superior bound for the domain wall 
width, corresponding to the case of single mode approxima- 
tion. In figures l6l and [TT1 this value was indicated by a hori- 
zontal line. We can see in each case how close profiles are of 
a sine function when the transition takes place. 



B. Model 2 

The mean field phase diagram for model 2 is shown in Fig. 
[7] An important difference with respect to the correspond- 
ing diagram of model 1 is the presence of a uniform phase 
for small T and a values. This diagram has some similarity 
with phase diagrams of anisotropic models with short range 
competing interactions, like the ANNNI model 2 ^. Neverthe- 
less, at variance with the ANNNI model, in the present case 
there is no direct transition from the uniform to the disordered 
phases, at least up to very small values of a, where the pre- 
cision of our algorithm breaks down. Instead, a sequence of 
two phase transitions with increasing temperature is observed. 
The modulated-uniform transition line characterizes a discon- 
tinuous transition at which the modulation length diverges, as 
can be seen in figure [8] Note that in the modulated region the 
free energy shows a single minimum at a finite modulation 
length A, but as the temperature is lowered a second mini- 
mum at the far right of the figure appears. Numerically, this 
minimum is compatible with an infinite value of A, typical of 
a homogeneous solution. 
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FIG. 6. (Color on-line) Temperature dependence of the domain wall 
width (£) normalized to the modulation length (A) for model 1 with 
a — 1. The superior full line (blue on-line) is the mean field solution 
and the lower one the solution with fluctuations included. The limit 
value of this magnitude for a perfect sine profile (single mode) is 
represented by the horizontal line. Dots represent the results of the 
numerical solution. 




a 



FIG. 7. (Color on-line) Mean field phase diagram for model 2. The 
dashed-dotted lines represent equal modulation length (A) curves . 
The full superior line (blue on-line) corresponds to a second order 
transition between a modulated and a disordered phase. The full 
lower line (red on-line) defines a transition between the modulated 
phase and a uniform phase. Dots represents the results of the numer- 
ical solution. 
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FIG. 8. (Color on-line) Free energy of model 2 as a function of mod- 
ulation length for different temperatures. Two top curves (blue on- 
line) show a single minimum at A slightly larger than 1 . In the two 
curves at the bottom (green on-line) a second minimum to the right of 
the figure appears, signaling a discontinuous transition to a uniform 
state . 
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FIG. 9. (Color on-line) Mean field phase diagram for model 2 in- 
cluding fluctuations. The dashed-dotted lines represent equal modu- 
lation length (A) curves . The top line (blue on-line) corresponds to 
a second order transition to a disordered phase. Below the full bot- 
tom line (red on-line) a uniform phase has less free energy than the 
modulated one. Note the bending of the constant A curves, meaning 
an increasing modulation length with a at the transition temperature. 
Dots represent the results of the numerical solutions. 



Figure [9] shows the phase diagram of model 2 when fluc- 
tuations to the mean field solution are included. As was the 
case for model 1, a first effect of fluctuations is to depress the 
transition line between the modulated and disordered lines. 

The curves of constant A for model 2 in the mean field ap- 
proximation are shown in Figure [7] They tend smoothly to the 
critical line, as in model 1 in such a way that higher curves 
corresponds to lower values of A. The modulation length at 
the transition is a constant independent of a, corresponding to 
that of the minimum of the fluctuation spectrum. On the other 



hand when fluctuations are included the shape of the equal A 
curves changes, as shown in Figure|9] While upper curves still 
corresponds to lower values of the modulation length, in this 
case the modulation length at the transition is an increasing 
function of a, at variance with the behavior of model 1 , where 
the modulation length at the transition is a decreasing function 
of the parameter a. 

In figure[T0]the behavior of the modulation length with tem- 
perature is shown for model 2. Note that the values of A stay 
very near Ao, the value at the transition, meaning that in the 
present model the modulation length has a weak dependence 
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FIG. 10. (Color on-line) Temperature dependence of modulation 
length (A) for model 2 with a — 0.5. Superior and inferior dots 
correspond to the results without fluctuations and when fluctuations 
are included, respectively. The full lines represent the best fits of 
numerical results by the proposed functionality d!7t . The lower line 
(red on-line) is interrupted because of the transition to the disordered 
phase.. 
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FIG. 1 1 . (Color on-line) Temperature dependence of the domain wall 
width (£) normalized to the modulation length (A) for a — 0.5 in 
model 2. The upper line (blue on-line) is the mean field solution and 
the lower one (red on-line) the solution with fluctuations included. 
The limit value of this magnitude for a perfect sine profile (single 
mode) is represented by the horizontal line. Dots represent the results 
of the numerical solution. 



with temperature. This is an important difference with respect 
to the behavior of model 1. Nevertheless, as shown by the 
solid lines in the same figure, a fit to the expression (fTTI i still 
works very well. 

The behavior of the domain wall width for model 2 is shown 
in figure [TT] Comparing with figure [6] for model 1, we note 
that in the present model the profile is always very close to a 
sine, even at low temperatures. Besides that, fluctuations seem 
to be unimportant for the behavior of the domain wall width, 
at variance with the behavior observed in model 1 . Neverthe- 
less, the amplitude of the order parameter behaves similarly to 
that of model 1, staying practically at saturation until the tran- 
sition temperature is reached, in the case when fluctuations are 
included. The fact that for model 2 the profiles stay very near 
the single mode case can be understood due to the presence 
of the quartic contribution in the fluctuation spectrum which 
implies that the energy required to excite higher harmonics at 
low temperatures is larger than in the quadratic case of model 
1. 



IV. ANOMALOUS BEHAVIOR OF STRIPE 
COMPRESSIBILITY 

A measure of a stripe compressibility in modulated systems 
can be obtained through the compression or Young modulus 
B(T). This magnitude is related with the energy cost associ- 
ated with deviations from the equilibrium modulation length 
(A(T)) and it is defined as±2: 



B(T) = A(T) 



2 5 2 F(A,T) 



8\ 2 



(20) 



\(T) 



We calculated the compression modulus directly from the free 
energy F(X, T). In a recent work, Portmann et al have pro- 
posed some scaling relations for modulated matter. A scaling 



relation for the compression modulus was obtained in terms of 
the equilibrium modulation length (A(T)) and the modulation 
amplitude (M(T)): 



B = cM(T) 2 A(T) A , 



(21) 



where c is a constant and A is an exponent related with the 
system dimensionality and the microscopic nature of interac- 
tions. One of the main predictions of the proposed scaling re- 
lations is an anomalous behavior of the Young modulus when 
the exponent A is negative. This is the case, e.g. in a dipolar 
frustrated ferromagnet in two dimensions, which has an effec- 
tive free energy of the type of model 1 in this work. More 
interestingly, Portmann and collaborators find the same form 
of scaling for the behavior of the critical magnetic field as 
a function of the modulation length^. This implies a reen- 
trant or inverse transition in the external field vs temperature 
phase diagram at low temperatures for all systems which have 
a negative value of the exponent A. This interesting result 
is supported by experimental measures in ultrathin films of 
Fe/Cu(001) from the same groupii, and it has been reported 
in a Ginzburg-Landau model for the same system^. 

In FiguresQ~2]and[T3]we show the results of a direct calcula- 
tion of the compression modulus (full curves) for both model 
1 and model 2 via equation (|20l . The typical behavior of 
B(T) in a normal system is a monotonically decreasing func- 
tion of temperature, growing temperature leading to weaker 
bonds between particles and a weaker elastic response. How- 
ever we can see in Figures[T2]and[T3]a regime for low enough 
temperatures where the compression modulus is an increas- 
ing function of temperature. Furthermore, the inclusion of 
fluctuations in both models do not destroy this behavior but 
strengthen the anomaly even more, which suggests a robust 
physics behind the numerical calculations. It is worth to men- 
tion that for very low temperatures the domain walls for model 
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FIG. 12. (Color on-line) Continuous curves corresponds to the calcu- 
lated Young modulus with and without fluctuations for model 1 . The 
dots represent the best fit using the scaling relation equation i2\\ . 



FIG. 13. (Color on-line) Continuous curves corresponds to the calcu- 
lated Young modulus with and without fluctuations for model 2. The 
dots represent the best fit using the scaling relation equation d2U . 



1 become very sharp and as result numerical instabilities ap- 
pear which can be seen in Figure Q~2] 

We also tested the proposed scaling relation by fitting our 
analytical results with expression (ETV Results are shown 
with dots in FiguresQ~2]and[T3] In all cases the proposed scal- 
ing law fits very well our numerical results. As part of our 
checking of the scaling predictions we let the initially fixed 
(quadratic) exponent of M (T) in equation d2lT i to vary. The 
results for this exponent as well as for A in model 1 were 
very close to those predicted ini2 for the 2d dipolar frustrated 
ferromagnet. On the other hand, for model 2 important devi- 
ations where found. Nevertheless, in this case deviations for 
the values of exponents are expected because, strictly speak- 
ing, model 2 does not correspond to the kind of systems con- 
sidered by Portmann and collaborators. 

It is interesting that the scaling relations seem to be valid 
for a wider class of models. Future work should address the 
extension of the scaling relations to more general microscopic 
interactions as well as the origin of the anomalous behavior of 
the Young modulus in stripe forming systems. 

V. CONCLUSIONS 

We have presented several new results regarding two well 
known effective models for stripe forming systems in two di- 
mensions. The models can be thought as the continuous limits 
of some microscopic models with competing isotropic inter- 
actions, or alternatively, they can be defined directly through 
the form of the high temperature structure factor. A com- 
mon feature of both models is that the structure factor has its 
maximum in a ring-like region in reciprocal space, which has 
strong consequences for the low temperature physics. Nev- 
ertheless, contrary to what is usually assumed in the litera- 
ture, important differences were found in their properties and 
phase diagrams that lie in the specific wave vector depen- 
dende of the high-T structure factor, with Model 1 defined by 
a non-analytic dependence on the modulus of wave vector and 



Model 2 defined by an analytic form. Model 1 is a good small 
k representation of some systems with long range interactions, 
like the dipolar frustrated ferromagnet and stripe forming sys- 
tems with elastic interactions with the subtrate, while Model 
2 can be obtained as an effective model for competing finite 
range isotropic interactions, which leads to an analytic expres- 
sion for the k dependence of the structure factor. 

We have computed the full phase diagrams of both mod- 
els in two approximations: mean field, and mean field plus 
fluctuations. Relevant equations were solved numerically. In 
both cases it is observed that fluctuations have a strong ef- 
fect on the phase diagrams and other observables. With re- 
gard to phase diagrams, the main effect of fluctuations is to 
change the order of the isotropic to modulated phase transition 
from continuous in the mean field approximation to a strongly 
discontinuous transition when fluctuations are taken into ac- 
count. Besides the isotropic-modulated transition, Model 2 
shows another discontinuous phase transition from the modu- 
lated to a uniform (infinite modulation length) phase below a 
critical value of the parameter a, which characterizes the cur- 
vature of the fluctuation spectrum at the minimum. 

In both models, the modulation length A is a smooth func- 
tion of temperature, characterizing incommensurate modula- 
tions. The temperature dependence of A is strong for Model 
1 and weak for Model 2. While in the mean field approxima- 
tion the modulation length attains its mimimum value at the 
transition, which corresponds to the wave vector at the mini- 
mum of the structure factor, when fluctuations are taken into 
account A can take values less than that corresponding to the 
minimum of the structure factor. Another quantity of interest 
is the width of the domain walls £. We have introduced an op- 
erational definition of the width, taking as a reference value 
that corresponds to a single mode, sine-like solution, valid 
near the transition in the mean field approximation. As was 
the case with the modulation length, £ has a stronger depen- 
dence with temperature for Model 1 . It is also important the 
fact that £ is smaller when fluctuations are included, meaning 
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sharper domain walls, in line with the sharp phase transitions 
obtained in this case, where the amplitude of the modulations 
stay near saturation already below the transition. 

Finally we have shown results for a response function, the 
compression modulus of stripes B(T). This quantity shows 
an anomaly in its behavior, when compared with normal sys- 
tems, having a maximum at intermediate temperatures. The 
anomalous behavior follows a recent prediction from a scal- 
ing hypothesis for modulated systems, which was obtained in 
the context of models with some particular long range inter- 
actions. We have verified its validity for both Model 1 and 
Model 2. Interestingly, although Model 2 does not belong to 
the kind of models for which the scaling hypothesis was pre- 
dicted, it shows the same kind of scaling, with a particular 
scaling exponent. The scaling exponent for Model 1 corre- 
sponds, to a good numerical approximation, to the value pre- 
dicted for a model with dipolar competing interactions. 

In summary, we have computed some relevant quantities 
for models with modulated phases in two dimensions, com- 
plementing old results and clarifiying some differences be- 



tween two models usually presented in the literature as giv- 
ing essentially the same physics. We have shown that, de- 
spite the obvious similarities, both models have important dif- 
ferences. An interesting outcome of our results is the evi- 
dence of anomalous behavior in a response function in both 
models considered, in line with recent results in related mod- 
els and experiments in ultrathin magnetic films. Future work 
should address the mechanisms that originate the anomalous 
behavior and the generality of this for other systems. Last 
but not least, the straightforward numerical implementation 
of the coarse-grained models defined in this work allowed us 
to access the whole temperature range of the phase diagrams, 
complementing already known results, usually limited to tem- 
peratures near the phase transition lines. 
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